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Abstract — We summarise the mathematical foundation of the 
holographic method of measuring the reflector profile of an 
antenna or radio telescope. In particular, we treat the case, where 
the signal source is located at a finite distance from the antenna 
under test, necessitating the inclusion of the so-called Fresnel 
field terms in the radiation integrals. We assume a "full phase" 
system with reference receiver to provide the reference phase. 
We describe in some detail the hardware and software imple- 
mentation of the system used for the holographic measurement 
of the 12m ALMA prototype submillimeter antennas. We include 
a description of the practicalities of a measurement and surface 
setting. The results for both the VertexRSI and AEC (Alcatel- 
EIE-Consortium) prototype ALMA antennas are presented. 

Index Terms — ALMA, Antenna measurements, radio hologra- 
phy, millimeter antenna, near-field, radio telescope. 



L Introduction 

LARGE reflector antennas, as those used in radio as- 
tronomy and deep- space communication, generaUy are 
composed of a set of surface panels, supported on three or 
more points by a support structure, often caUed the backup 
structure. After assembly of the reflector it is necessary to 
accurately locate the panels onto the prescribed paraboloidal 
surface in order to obtain the maximum antenna gain. The fact 
that some antennas have a "shaped"contour is irrelevant for the 
purpose of our discussion. We are concerned W\\h describing a 
method which allows us to derive the position of the individual 
panels in space and compute the necessary adjustments of 
their support points to obtain a continuous surface of a certain 
prescribed shape. 

The analysis by Ruze [1] of the influence of random errors 
in the reflector contour on the antenna gain indicates that the 
RMS error should be less than about one-sixteenth of the 
wavelength for acceptable performance. Under the assumption 
that the errors are small compared to a wavelength, randomly 
distributed with RMS value e, have a correlation length c 
which is much larger than the wavelength A, and much smaller 
than the reflector diameter D, the relative decrease in aperture 
efficiency (or gain) can be expressed by the simple formula 



exp 




(1) 



where r]Ao is the aperture efficiency of the perfect reflector. 
An error e of A/40 is required to limit the gain loss to 10 
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percent; with an error of A/16 the gain is decreased to about 
half of the maximum achievable. 

The setting of the reflector panels at accuracies better than 
100 fim has required the development of measuring methods 
of hitherto unsurpassed accuracy. It should be noted that these 
measurements need to be done "in the field", which in the 
case of millimeter radio telescopes generally means the hostile 
environment of a high mountain site. One versatile, and by 
now widely used method is normally called "radio hologra- 
phy". The method makes use of a well-known relationship 
in antenna theory: the far-field radiation pattern of a reflector 
antenna is the Fourier Transformation of the field distribution 
in the aperture plane of the antenna. Note that this relationship 
applies to the amplitude/phase distributions, not to the power 
pattern. Thus, if we can measure the radiation pattern, in 
amplitude and phase, we can derive by Fourier Transformation 
the amplitude and phase distribution in the antenna aperture 
plane with an acceptable spatial resolution. Bennett et al [2] 
presented a sufficiently detailed analysis of this method to 
draw the attention of radio astronomers. Scott & Ryle [3] used 
the new Cambridge 5 km array to measure the shape of four of 
the eight antennas, using a celestial radio point source and the 
remaining antennas to provide the reference signal. Simulation 
algorithms were developed by Rahmat-Samii [4] and others, 
adding to the practicability of the method. Using the giant 
water vapour maser at 22 GHz in Orion as a source Morris 
et al [5] achieved a measurement accuracy of 30 /im and were 
able to set the surface of the IRAM 30-m nuUimeter telescope 
to an accuracy of better than 100 //m. 

Artificial satellites, radiating a beacon signal at a fixed 
frequency have also been used as farfield (Rf = ^y-) 
signal sources. Extensive use has been made of synchronous 
communication satellites in the 11 GHz band [6], [7]. These 
transmitters of course do not provide the range of elevation 
angles accessible with cosmic sources. Some satellites, notably 
the LES (Lincoln Experimental Satellite) 8 and 9, have been 
used for radio holography of millimeter telescopes [8]. They 
provided a signal at the high frequency of 37 GHz and with 
their geo- synchronous orbit moved over some 60 degrees in 
elevation angle. Unfortunately, both satellites are no longer 
available. Lacking a sufficiently strong source in the farfield, 
we have to take recourse to using an earth-bound transnutter. 
In practice these will be located at a distance of several 
hundreds of meters to a few kilometers and be at an elevation 
angle of less than 10 degrees. Clearly, these are in the nearfield 
of the antenna, requiring significant corrections to the received 
signals. In particular, the phase front of the incoming waves 
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will not be plane and it contains higher order terms in the 
radial coordinate of the antenna aperture. These must be 
corrected before the Fourier transformation can be applied. 
We treat these corrections in detail in this paper. 

Successful measurements on short ranges have been re- 
ported for the University of Texas millimeter telescope [9], 
the IRAM 30-m telecope [10], the JCMT [11] and the ASTE 
antenna of NAOJ [12]. 

ALMA (Atacama Large Millimeter Array) is a new large 
aperture synthesis array for submillimeter astronomy consist- 
ing of 50 high accuracy antennas of 12 m diameter. The 
instrument is under construction at 5000 m altitude in the 
Atacama desert of northern Chile. ALMA is a collaboration of 
North America and Europe with participation of Japan. Two 
prototype antennas were procured and erected at the site of 
the Very Large Array of NRAO in New Mexico. The results 
of an extensive evaluation program of these antennas has 
been presented by Mangum et ah [13]. The reflector surface 
accuracy was specified at 20-25 /im, requiring a measurement 
method with an accuracy of 10 /im or better. This was achieved 
with a near-field holographic system using a transmitter at a 
wavelength of 3 mm and at a distance of only 315 m from the 
antennas at an elevation angle of 9 degrees. Here we describe 
these measurements in some detail. 

IL The Mathematics of Radio Holography 

The reciprocity theorem describes the equivalency between 
the characteristics of a transmitting and receiving antenna. 
Thus both concepts will be used in the following treatment 
depending on the specific aspect under description. We shall 
not repeat here the fundamental analysis which leads from 
Maxwell's equations to the "physical optics" representation of 
the characteristics of the reflector antenna {e.g. [14], [15]). 
Our discussion below essentially follows the treatment by 
Silver [14]. The basic expression, linking the radiation function 
f{x^y^z) at a point P in space with the field distribution 
7^) over the aperture plane of the antenna, is written as 
(see Fig. [T] for the geometry) 



f{x,y,z) = J 



47rr 



-ikr 



ik - 



iz ' r-i^ ikiz • s 



(2) 



where the integration is extended over the aperture area, 
k = 2ti/\ and the unit vectors are as indicated in Fig. [T] (with 
s the propagation vector of the wave field in the aperture). 
This relation assumes that the aperture is large in units of 
the wavelength. This general expression can be simplified 
depending on the distance of the field point P from the aperture 
plane. We discern the so-called far-field region (Fraunhofer 
diffraction), near-field region (Fresnel diffraction) and the 
evanescent-wave zone up to a few wavelengths from the 
reflector. In the evanescent- wave zone, which does not concern 
us here, no approximations are allowed. 

If the field point P is sufficiently far away from the aper- 
ture, the following simplifications can be introduced in the 
evaluation of Eq. [2j 



PCs.y.z) 




Fig. I. Geometry of the aperture integration method for finite distance to 
the field point P. 



1) The term ^ is ignored with respect to k in the bracketed 
term. 

2) The term ^ outside the brackets is replaced by the 
reciprocal distance from the aperture center to the 
field point P. 

3) The term iz - Ti can be approximated by iz - Ri = 
cosO with i^i the unit vector from the origin to the 
field point. 

4) the term iz - s represents a deviation from uniform 
phase over the aperture. If these are small, this term 
can be assumed to be equal to one over the aperture. 

With these approximations Eq. |2] is simplified to 



f{x,y,z) 



2XR 



J F(^, T]) [cos (9 + 1] e'^'^d^df]. (3) 

For the distance r from any point in the aperture to the field 
point P we have (see Fig. [T]) 



{{x-0^ + {y-vf + z'} 



0.5 



(4) 



Writing the coordinates of the field point P(x,y,z) in spherical 
coordinates, we obtain 



X = R sin cos (f) = Ru, 
y — R sin sin (j) = Rv , 



z = RcosO = R^(l-v? -v'^), 



(5) 



where we have also introduced the direction cosines of the 
field point 

P(ii, v) = (sin cos (/), sin 9 sin <p) . 
Thus, Eq. |4] can be written as 



r = {{Ru - ^ {Rv - r]f ^ R\l - - v^)}""'^ 



R<l-2 



R 



R^ J 



(6) 
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The series expansion of Eq. [6] yields 



R — {u£_ + VT]) 

« + vr]f 
2R 



2R SR^ 



(7) 



A. The Far-Field Approximation (Fraunhofer Region) 

In the far-field situation, where R tends to infinity and 
Ri and ri are essentially parallel, the variation of r in the 
exponent of Eq. [3] can be reduced to the linear part of Eq. [7] 



r = R — {u^ + VT]) 



(8) 



B. The Near- Field Approximation (Fresnel Region) 

In the near-field region, which corresponds to the Fresnel 
region in optical diffraction (e.g. [17]) most of the simplifi- 
cations leading to Eq. [3] can still be used, as long as r is at 
least several aperture diameters large. However, the variation 
in r over the aperture must now include higher-order terms in 
Eq. [7] and be maintained in the exponent (phase) term of the 
integral. Normally, for the Fresnel region analysis, the series 
is stopped after the quadratic term, which preserves the first 
three terms of the series in Eq. [7] Thus the near-field (Fresnel 
region) expression can be written in the form of the following 
radiation integral, which is the well-known Fresnel diffraction 
integral in two coordinates: 



Additionally, considering that for a high gain antenna the 
angular region of interest is confined to small values of 0, 
we can write Eq. [3]with cos^ = 1, which is valid to 0.1% for 
angles up to 3 degrees off-axis. The radiation integral (Eq. [3]), 
then becomes 



^ikR 



A R 



exp < ik 



-{U^ + VT]) + 



2R 



■d^df]. (11) 
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ikR 



R 



F{£^^r]) exj[){—ik{^u -\- r]v)}d^dr] (9) 



where the integration is performed over the aperture A. We 
see from Eq. [9] that there exists a Fourier Transformation 
relationship between f{u^v) and F(^,7^). Ignoring the term 
J, the inverse Fourier transformation can thus be written as 



In the application of holography in the near-field we want 
to derive the complex aperture field distribution from the 
measured near-field pattern. Thus the inverse Fourier Trans- 
formation of Eq. [3] will be our point of departure, where Eq. [6] 
is the expression for the finite distance from a point in the 
aperture to the point where the signal source is located. Thus 
we have the inverse of Eq. [3] 



1 AkR r ^ f 

Fi^^v) = J f{u^v)eMik{u^^vr])}dudv, (10) F{^.v) = J 



47r R 

where the integration in principle has to be performed over 
a closed surface, surrounding the aperture. Thus a knowledge 
of the entire far-field pattern f{u^v) both in amplitude and in 
phase provides a description of the complex field distribution 
F(^, 7^) over the aperture plane of the antenna, also in ampli- 
tude and phase. This forms the basis of the so-called radio 
holographic measurement of the shape of a reflector antenna. 
Deviations from a uniform phase function over the aperture 
are thereby linked to local errors in the prescribed contour of 
the reflector surface. 

It is interesting to note that Silver devotes a lengthy dis- 
cussion to this Fourier Transform relationship ([14], Ch. 6.3), 
but concludes that the practical application is limited by the 
fact that the far-field pattern is only prescribed in power. 
Thus the phase function of f{u^v) would be arbitrary and 
the aperture distribution cannot be uniquely determined. It 
was Jennison [16] who mentioned the same relation and its 
practical usefulness, pointing out that the amplitude and phase 
can both be measured with an interferometer. When Silver 
wrote his text in the mid- 1940s, radio interferometry had not 
yet been developed. 

In most cases it will be impossible, or in any case impracti- 
cal, to measure the far-field pattern over the entire sphere. The 
Nyquist sampling theorem shows however that a measurement 
of the pattern out to an angle 9 = n 6^ from the beam 
axis yields the aperture distribution with a spatial resolution 
of S = ^, where 9^ ~ ;^ is the half-power beam- width, D 
is the aperture diameter, and A the wavelength. 



f{u^ v) ex.p{—ikr)dudv . 



(12) 



where R is the distance from the antenna aperture center to 
the holography signal source. We maintain all terms of Eq. [7] 
in order to make an estimate of the error with respect to the 
usual Fresnel approximation. 
We rewrite Eq. [7]as follows: 



r ^ R - [u^^vT]) ^ Spi (^, 7?) + e 



(13) 



where we define the terms, which are independent of the 
integration variables, as the variable Sp^: 



2\2 



2R 



(14) 



while the other terms in higher powers of (w, v) are collected 
under the variable e. 



Substitution of Eq. |7] into Eq. [12] yields 



(15) 



-ikR 



■exp{-iA:(5pi(^,7?)} 



A R 

f{u^ v) ex.p{ik{u^ + vr])}e~^^^dudv, (16) 
The terms in e "modify" the direct Fourier transformation of 



Eq. 16 





Fig. 2. The geometry of axial (upper half) and lateral (lower half) 
displacement of the feed from the focus of a parabola. 



1) The 6pi Term: The first path-length term 5pi causes a 
rapidly varying phase variation over the aperture, which can 
be compensated to a large degree by an axial displacement 
of the feed. In Fig. [2] we illustrate the geometry of both a 
lateral 6fi and an axial 6f^ displacement of the feed from the 
primary focus. We need to calculate the difference in path- 
length {p' — p) as a function of the radial position r of the 
point P. The residual path-length error at P with respect to the 
ray to the vertex V will then be given by 

= [p' -P)p-{p' -P)y^ (17) 

where = + Applying Pythagoras' law to the upper 
triangle PQF' and using the defining relation for the parabola 
= where z = VQ, we find for the path-length varia- 
tion due to an axial defocus 5f away from the reflector 




0.5 



(18) 



We want to minimise the sum of the two terms {5pi + 5p2) 



(Eqs. [14] and [18]) by choosing the appropriate value of 3f. 
Because of the (^,7/)-dependence there will be a residual path- 
length error, which we must apply to the result of the Fourier 
Transformation. Fig. [3] shows the residual path-length error 
{5pi + 5p2) for several choices of 5f for the geometry of the 
ALMA antennas and the actual distance to the holography 
transmitter. A value of 96-98 mm limits the error to ±3 mm 
over the aperture for R = 315 m. This remaining error must be 
introduced in the mathematical analysis of the data according 



to the curve shown in Fig. [3] This is a correction to the aperture 
phase distribution, obtained after the Fourier Transformation 
of the measured beam pattern. 



2) The e Term: The higher order terms in Eq. 15 (e), 
containing the integration variables iu^v), constitute a small 
path-length error which adds a phase term to the integral in 
Eq. [16] of the form 



exp(— i/ce) 



1 — ike 



1 — ik lu 



^2 



2^ 2V' 

u —- — V —- — uv— 
2R 2R R 



} 



(19) 



It is seen that this correction involves the calculation of five 
additional integrals, which look like Fourier transformations, 
but aren't really bona fide Fourier transformations. When all 



the integrals of Eq. 19 are evaluated, it turns out that the 
contribution of e amounts to 2pm of path-length over most of 
the aperture, reaching a peak value of 10/im at the very edge 
of the aperture. This is illustrated in Fig. [4] In a high accuracy 
measurement, where the aim is to achieve a measurement 
accuracy of better than 10 /im, it is advisable to include this 
correction term. 

3) Dependence on Aperture Plane Reference: In the cor- 
rections for the finite distance to the transmitter (the near- 
field corrections), we define the aperture plane at a convenient 
location, normally halfway between the vertex and the edge of 
the reflector. We have taken the center of this aperture plane 
as the origin of the coordinate system. In most antennas there 
is a significant distance between this plane and the axes of 
rotation for the movement of the reflector (see Fig. [5]). From 
this figure we see that there is a "parallax" effect between the 
adopted direction cosines (u, v) and those given by the antenna 
scanning coordinates (u^ vO, given by the relations 
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Fig. 4. Non-Fresnel correction terms over the aperture of the 12 m diameter 
ALMA antenna. Horizontal and vertical axes units are meters with color 
showing surface error in jim shown at right. The departures from circular 
symmetry are due to the effects of the actual surface errors present in the 
map. 



-100 




Fig. 5. Illustration of the geometry of selected aperture plane and antenna 
rotation axis. T is the position of the transmitter while scanning the antenna. 



radlus[m] 



Fig. 6. Top: The near-field path-length correction as a function of radius 
assuming A = 3.1 m and R = 315 m. Bottom: The surface error resulting 
from the misregistration of the radial coordinate {i.e. ignoring the difference 
between direction cosines u and u'). 



significant, causing a surface error as a function of radius as 
shown in the lower part of Fig. [6} its RMS value is 18/im in our 
case, significant with respect to the required setting accuracy. 
Such a donut- shaped systematic deviation was indeed found 
in our original maps, once we applied the geometry correctly. 
It was of course treated properly in the final measurements 
and surface setting. 

4) Dependence on Focal Deviation: It is possible that 
during the measurement the receiver feed is not located in the 
optimum focal position. With reference to Fig. [2] and Eq. 18 
it can be shown that the path-length error caused by an axial 
defocus of 5z is given by 



u' (^1 + 



A' 
A' 



(20) 



where A is the distance between rotation axis and aperture 
plane. We use the scanning coordinates (u^ vO, read from the 
antenna encoders, to calculate the position of the points in the 



aperture plane. If (u, v) in the Fourier integral (Eq. [16]) are 
approximated by {u' .,v'), the scale of the aperture map will 
be overestimated by the factor (l + The result of this is 
that the near-field correction for each pixel in the map is not 
evaluated at the correct radius. This causes a path-length error 
proportional to the derivative of the near-field correction with 
respect to the radial coordinate. Fig. |6] illustrates the magnitude 
of this effect for the case of our geometry, where A ^ 3.1m, 
i.e. about 1% of the distance R to the transmitter. The error is 



5p^ = 5z{l 



1 



4/^ 



51 
f 



(21) 



4/2 I V" 4/2 

while a transverse (lateral) offset by an amount (5x in the ^- 
plane (Fig. [2j lower half) will cause a path-length variation 
of 



/ 1 1 



4/2 



(22) 

In the reduction process of the holography data, these terms 
are found by a fit of the measured beam map. The final map 
of surface deviations is then referred to a position of the feed 
in the fitted "out-of-focus" location. 
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Fig. 7. The VertexRSI (left) and AEC (right) ALMA prototype antennas. 

III. PRACTICAL Application of the Holography 
Measurements 

A. Task 

We now describe the way in which a holography measure- 
ment has been executed on the ALMA prototype antennas 
(Fig. [7]). The specification requires the antennas to have a 
surface accuracy of 25/im RMS for the AEC antenna (with a 
goal of 20/im) and 20/im for the VertexRSI antenna. ALMA 
assumed the task to demonstrate this with the aid of a 
holography system at 3 mm wavelength after delivery of the 
antennas by the contractors with a surface accuracy of not 
worse than 100/im RMS. This initial setting was performed 
by VertexRSI with digital photogrammetry and by AEC with 
the aid of a Leica "total station" laser-tracker (essentially a 
theodolite with integrated distance measurement instrument 
and all-electronic readout). 

The holography system was designed to provide a measure- 
ment repeatability of 10/im, which would suffice to demon- 
strate the realism in the obtained overall surface accuracy. It 
should be noted that in the current setup the holography system 
provides a surface map at one elevation only. No information 
on the gravitational deformation of the antenna with varying 
elevation angle can be obtained. 

B. Equipment and Measurement Program 

• The signal source for the holography measurements is 
a monochromatic transmitter at a frequency of 78.92 or 
104.02 GHz, located on a 50 m high tower at a distance 
of 315 and 302 m from the VertexRSI and AEC antenna, 
respectively. The elevation angle is approximately 9 de- 
grees. 

• The receiver is a full-phase double-receiver, located in 
the apex region behind the primary focus of the main 
antenna. The reference signal is received by a wide beam 
horn pointing along the boresight towards the transmitter. 

• Amplitude and phase maps of the antenna beam were 
obtained by raster scanning. The Nyquist sampling theo- 
rem provides the link between the angular size of the 
observed map and the required spatial resolution over 
the aperture. If we want to obtain n independent samples 



over the diameter of the aperture, we need to extend the 
map to an angle of n times the half-power beamwidth 
off axis. We chose a map size to obtain about 0.15 m 
spatial resolution after Fourier transformation of the map. 
A typical measurement then takes about one hour of time. 

• From the phase distribution, which is a representation of 
the misaUgnment of the 264 (VertexRSI) or 120 (AEC) 
panels constituting the reflector, the necessary adjust- 
ments of the 5 support points per panel were derived. 
These were then applied by hand with a simple tool to 
improve the accuracy of the reflector surface. 

The algorithms and software used for the data analysis 
and derivation of the panel adjustments have been applied 
successfully at the telescopes of IRAM. The necessary correc- 
tions for the finite distance to the transmitter (the "near-field" 
corrections) in our case were derived and checked against 
similar corrections applied by others, e.g. for the JCMT [11]. 

The equipment has been designed to provide sufficient 
signal-to-noise ratio to render the error due to noise insignif- 
icant. The greatest risk in this type of measurement lies in 
undetected or poorly corrected systematic errors. 

• An accurate knowledge of the amplitude and phase 
function of the feedhorn, illuminating the reflector, is 
essential, because errors in these are fully transferred to 
the aperture phase map and hence to the surface profile. 

• Multiple reflections from the ground or structures form a 
possible source of errors in this type of work. We care- 
fully covered all areas of potentially harmful reflections 
with absorbing material. In some controlled experiments 
we could not demonstrate the existence of reflections. 

• The dynamic range of the receiver must be sufficient to 
accommodate the strong signal on the peak of the beam 
and the very weak signals towards the edge of the scan. 
There might have been some saturation on some of the 
measurements. We discuss this in more detail below. 

• The effect of the finite distance of the transmitter can 
be removed to a large extent (but not completely) by an 
axial shift in the position of the feed. An error in the 
distance to the transmitter thus can be corrected in the 
data analysis by a small adjustment of the feed position. 
The remaining phase error can be accurately calculated 
and applied to the data. 

C. Holography System Hardware 

The hardware specifications and requirements are sum- 
marised in Tables [I| and [Il| In the following we briefly describe 
the hardware components that comprise the holographic mea- 
surement system. 

1 ) Front-End: The front-end (see Fig. [8| is enclosed in a 
small, temperature controlled box with a diameter of about 30 
cm and a length of 50 cm. It fits inside the "apex structure" 
behind the primary focus of the VertexRSI antenna. The AEC 
antenna does not provide such a wide space and the receiver is 
bolted to the outside flange of the apex structure with a long 
piece of waveguide bringing the signal feed in focus. Both 
the signal- and reference-receiver are housed "back- to-back" 
in this box. This provides a compact system in which the LO 
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TABLE II 

Holography Hardware Specifications 



Frequencies 

Frequency Stability 

Receiver Bandwidth 

Receiver Tunability 

Transmitter Antenna Gain 

Transmitter EIRP (P) 

Transmitter Power to Antenna 

Transmitter Antenna Beam- Width @ — 3dB 

Reference Antenna Beam- Width @ — 3dB 

Main Feed Beam- Width @ -3dB 

System Temperature 

Reference Feed Power Received (Pr) 

On-Boresight Signal (Mq) 

On-Boresight Noise ((jq) 

Off-Boresight Noise (Pr Term) 

Average map noise for complex correlator {(Tav) 



78.92 and 104.02 GHz 

< ±5 Hz/day 

10 kHz 

130 MHz 

33dB 

> 20mW 

> lOnW 

4.6 deg (twice antenna angle at xmtr) 

4.6 deg (twice scan range) 
128 deg (-3dB edge taper) 
3200 K 

1.7 X 10-9p 
4.2 X lO-'^P 

(1.2 X 10-22l^(p))i 
(2.1 X 10-2'^ H^(P))^ 
(2.2 X lO-^^l^(P))^ 



TABLE I 

HOLOGRAPHY HARDWARE REQUIREMENTS 



Measurement Error 


< 


10/im 


Phase Accuracy 


< 


0.3 deg (2.5/im @ 3mm) RMS 


Amplitude Accuracy 


< 


1% 


Dynamic Range 


> 


43dB 


Signal-to-Noise Ratio (SNR) 


> 


40dB 


Channel-to-Channel Isolation 


> 


lOOdB 


Data Rate 




80 samples/sec (12 msec) 




Fig. 8. Holography system hardware. Left: Signal feed side of the front-end. 
Middle: Reference feed side of the front-end. Right: transmitter on top of 
the tower pointing at the VertexRSI prototype antenna in the foreground. The 
other antennas are part of the NRAO Very Large Array (VLA). 

signals can easily be made equal in length, greatly contributing 
to the phase stability of the system. Broadband mixers at 
ambient temperature convert the received signal frequency 
to a baseband of 10 kHz width. The system is designed 
with two frequencies at 78.9 and 104.02 GHz. Making the 
measurement at two different frequencies can be helpful in 
discerning systematic effects in the resulting maps, for instance 
caused by multiple reflections. The receiver is also tunable 
around each of these frequencies by 130 MHz for similar 
reasons. The signal horn is a conical, grooved cylindrical 
waveguide horn, while the reference horn is of similar design 
and equipped with a lens to provide a reference beam with a 
beam- width of 4.6 degrees at the half-power points. 

As is clear from the theoretical treatment above, it is 
imperative that we know the amplitude and phase function of 
both the reference and the signal feed as accurately as possible. 
The phase function must be subtracted from the measured 
aperture phase before connecting its phase variations to errors 
in the reflector profile. The feedhorns have been measured 
with great care on the indoor range at IRAM in Grenoble 
[18]. The results were compared with model calculations using 



an advanced electro-magnetic simulation package (the FDTD 
package of Microwave Studio from Computer Simulation 
Technology) and excellent agreement was found. The phase 
pattern of the feeds have an estimated error of less than one 
degree, while the amplitude taper at the edge of the reflector 
aperture is —6 dB. This is more than we would like (a free- 
space taper of 2.5 dB has to be added to the measured level). 
For a high signal to noise ratio in the outer part of the reflector 
an actual level of — 6 dB is preferred. For the measurement of 
the ALMA production antennas this feed should be replaced 
by one which provides such a taper. 

2 ) Back-End and Transmitter: The back-end of the receiver 
is essentially a digital signal processor (DSP) where the 
narrow-band signals are digitized and correlated. Both the 
"sine" and "cosine" part of the complex correlation function 
are obtained, which are then transformed to the amplitude and 
phase functions. 

The transmitter consists of a single photo-diode, directly 
coupled to a waveguide horn, which is fed through an optical 
fiber by two optical signals at different frequencies near a 
wavelength of ~ 1550 nm. The photo-diode provides a mixing 
signal at the difference of the two optical signals, tunable 
roughly from 78.7 to 79.0 GHz (low band) and 103.8 to 
104.2 GHz (high band), with an output power of about 10 
nW, leading to an EIRP of about 20 /iW. The transmitter is 
placed on top of a 50 m high tower at a distance of 300 to 325 
m from the aperture of three antennas at the site, resulting in 
a measurement elevation angle of about 9 degrees. 

IV. Holographic Data Acquisition 

To derive typical values for the various holography map 
parameters, we set the following boundary conditions: 

• The data rate is the canonical 12 msec per sample, which 
means about 80 samples per second. 

• The fine tuning feature of the holography receiver allows 
for the search for ground reflection. 

• A goal for the total time for one map is less than one 
hour. 

• The required aperature plane resolution is < 20 cm. 
This yields > 25 independent points per square meter 
of reflector surface. 
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• The 5 panel support points are on average some 0.4 to 0.6 
m apart for the VertexRSI and AEC panels, respectively. 
A twist in the panel with a similar scale length can be 
partially corrected. With a measurement spatial resolution 
of 15-20 cm this large scale twist can be fitted sufficiently 
well. 

• Oversample by a factor of at least 2 to minimize aliasing. 
Based on the equations listed in Appendix |lj with /i = 

1.13(6 + 2.5 dB taper), we obtain the typical holography map 
parameters of Table |IIl| 

V. Holographic Data Analysis 

A. Description 

Data analysis uses the CLIC data reduction software of the 
Plateau de Bure interferometer. The raw data, written by the 
on-line software in the ALM ATI-FITS data format [19], is 
converted to Plateau de Bure format using CLIC. The data 
are then calibrated and imaged using CLIC. The two main 
operations are: 

1) Calibrate data in amplitude and phase, based on bore- 
sight measurements at beginning and end of each map 
row, assuming gradual drift in amplitude and phase 
with time. This uses the standard amplitude and phase 
calibration commands in CLIC, which: 

a) Fit cubic spline functions of time to the observed 
amplitude and phase data on the boresight mea- 
surements. 

b) Subtract the phase spline function from the ob- 
served phase for the mapping scans. 

c) Divide the observed amplitude for the mapping 
scans by the amplitude spline function. 

2) Compute the aperture map and fit panel displacements 
and deformations. The data processing steps for compu- 
tation of the aperture maps are: 

a) Interpolate data to a regular grid in the antenna- 
based coordinate system. This grid matches the 
observed system of rows (same number and sepa- 
ration). This grid is further extended, by addition 
of zeroes, to a user-specified size, in order to get 
a finer interpolation of the output aperture map: 
64x64, 128x128, 256x256 and 512x512 sizes are 
available. 

b) FFT to aperture plane. This is replaced by a more 
complex transformation if one takes into account 
the first non-Fresnel terms, as described in ^ |II-B 



(Eq.[T6]). 

c) Compute phases in the aperture plane. 

d) Apply the geometrical phase correction: This is 



Eq. 14 plus Eq. 18 substituting p — + r]2 
as the radius in the aperture, / the focal length 
of the primary, and Sf by the distance between 
the holographic horn phase center and the antenna 



prime focus (see §II-B ). 
e) Correct for measured feed phase diagram. 

The measurement is described in the memo by 
Lazareff et al. [18]. 



f) Mask edges and blockage. 

g) Fit and remove 6 phase terms: Constant, 2 linear 
gradients in the horizontal and vertical directions, 3 
focus translations. These terms account for a phase 
offset, an antenna pointing error (constant during 
the measurement) and a small vector displacement 
of the holography horn relative to the nominal 
focus position (/ + df). The phase terms for the 
axial and transverse displacements of the focus are 
given in Eqs [21] and [22] respectively. Optionally one 
may keep fixed either the 5x and 5y coordinates 
or all three 5x, 6y, and Sz coordinates. 

3) Convert to normal displacement map. See Appendix |II| 
for details. 

4) Plot amplitude and phase maps. 

5) Fit panel displacements (optionally deformations) and 
screw adjustments. 

a) Each panel is assumed to be displaced (in the axial 
direction), tilted (around two orthogonal axes), 
and possibly deformed (deformation is a quadratic 
function of position offset relative to the panel 
center). As there are only five screws, only two 
deformation modes are allowed, we have thus five 



(23) 
(24) 
(25) 
(26) 
(27) 

where x and y are coordinates of a point on a panel 
(x and y in the plane tangent to the panel surface, 
with X axis in the radial direction). 

b) The five coefficients a, 6, c, e are independently 
fitted for each panel to the relevant part of the map 
for this panel. 

c) To take into account the effects of finite angular 
resolution, an iterative procedure is used where: 

i) The radiated beam is calculated using the fitted 
panel surfaces. 

ii) This beam is truncated to the size of the 
observed beam map, and subtracted from the 
observed beam. 

iii) An incremental surface map is calculated. 

iv) A new set of incremental panel displacements 
is calculated from this map. 

The procedure converges after a few iterations. The 
screw settings are output to a text file. 
6) These screw settings are applied to the panel adjusters 
to improve the surface accuracy of the reflector. The 
adjustments were done with a simple tool. Two people 
on a manlift approached the surface from the front, 
where the adjustment screws are located (see Fig.|9]). The 
time needed for an adjustment of the total of 1320/600 
adjusters was 8/7 hours for the VertexRSI/AEC pro- 
totype antennas, respectively. The entire procedure is 



displacement modes: 






Spi{x,y) 


= a 




Sp2{x,y) 


= bx 




^Psix.y) 


= cy 




Sp4{x,y) 


= d{x^ 




^Pb{x,y) 


= e{x^ 
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TABLE III 

Typical Holography Map Parameters 



Map Type 


5d 


fosr 


Oext 


sr 





Nrow 


foss 


tmap 




(cm) 




(deg) 


(arcsec) 


(arcsec/s) 






(hr) 


Standard 


20 


2.2 


1.64/1.24 


33/25 


300 


180 


20/15 


0.96/0.73 


Fine 


13 


2.2 


2.46/1.87 


33/25 


600 


270 


40/30 


1.08/0.82 


Course 


20 


1.4 


1.64/1.24 


53/40 


300 


112 


20/15 


0.61/0.46 



Assumes fi = 1.13 (6 + 2.5 dB taper), ly = 78.92/104.02 GHz, Ob = 74/56 arcsec, and /, 



1.3. 




Fig. 9. Panel adjustment of the VertexRSI prototype antenna. 

repeated until the required accuracy is achieved, or 
alternatively until the inherent measuement accuracy has 
been reached and no further improvement is surface 
accuracy is obtained. 



-100 -50 50 100 




Fig. 10. Final holography maps for the AEC (left) and VertexRSI (right) 
ALMA prototype antennas. Horizontal and vertical axes units are meters 
with color scale plotted as surface displacement in /im (top right). Note that 
the "x" (AEC) and "+" (VertexRSI) feed leg structure and the difficulties 
encountered with holographic measurements near these structures leads to the 
poor measurement results in these areas. 



VI. ALMA Prototype Antenna Holography 
Measurement Results 

The holographic measurement and setting of both anten- 
nas was performed immediately after the antennas became 
available for evaluation. During a period of about one year 
the antennas were subjected to a number of hard loads, like 
fast switching tests, drive system errors resulting in strong 
vibrations, and high-speed emergency stops collisions. Also, 
the influence of wind and diurnal temperature variations on 
the surface stability was a point of concern. It was therefore 
decided to close the evaluation program with a second holo- 
graphic measurement of the reflector surfaces. This was done 
in December 2004 to February 2005 during relatively good 
atmospheric conditions. It was during these measurements that 
we discovered that we had not correctly taken care of the 



correction explained above in ^ II-B.3 This correction was 
subsequently applied properly and the final surface maps are 



shown in Fig. [TO] Both antennas were set to a surface accuracy 
of 16-17 fim RMS. 

A. VertexRSI Antenna 

1) Overview: The antenna was delivered with a nominal 
surface error of 80 /im RMS, as determined from a pho- 
togrammetric measurement. Our first holography map showed 
an RMS of approximately 85 /im. A first setting of the 
surface resulted in an RMS of 64 /im. In four more steps of 
holographic measurement followed by adjustment the surface 
error decreased to 19 /im RMS. 

The sequence of surface error maps, along with the RMS 



and the error distribution is shown in Fig. 11 As allowed 



in the specification, we have applied a weighting over the 
aperture proportional to the illumination pattern of the feed. 
This essentially diminishes the influence of the surface errors 
in the outer areas of the reflector. The white areas in the surface 
error maps are the quadripod, optical pointing telescope, and 
a few bad panels, which could not be set accurately. All of 
these structures were left out of the calculation of the final 
overall RMS value. 

With increasing accuracy the presence of an artefact in 
the outer area of the aperture became apparent. There is a 
"wavy" structure in the outer section with a "period" too 
large to be inherent in the panels. Experiments with absorbing 
material showed that it was not caused by multiple reflections. 
The effect can be described by a DC-offset in the central 
point of the measured antenna map, i.e. some saturation 
on the point with the highest intensity. By adjusting this 
offset in the software, most of the artefact could be removed. 
This has been done with the final data. The additional set 
of follow-up holography measurements in December 2004 - 
February 2005 did not suffer from this signal saturation, and 
no artefact was observed in these follow-up maps. Checks 
of the holography hardware indeed suggest that the 2003 
holography measurements did experience a small amount of 
signal saturation. 

The best surface maps were obtained at night. During the 
spring 2003 period they consistently show an RMS of about 
20/im. Daytime maps tend to be somewhat worse; typical 
values of the RMS lie between 20 and 25/im. Part of this 
is certainly due to the atmosphere, even over the short path- 
length of 315 m. 
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28-apr-2003 
67.6 fim 



02-may-2003 
38.6 fxm 



07-may-2Q03 
28.3 fxrc] 




08-may-2J03 
21 .6 Mm 



15-may-2q03 
19.1 fxm 



10-feb-20Q5 
16.1 fxm 



Surface Error Map (/^m) 



-200-100 100 200 

Error Distrib. (%) vs. Error (/^m) 



Fig. 11. Sequence of surface error maps with intermediate panel setting. 
The surface contours are shown on the left side; the error distribution, plotted 
as percentage contribution versus surface error, on the right. The white cross 
and the small white areas represent the quadripod and a few faulty panels 
and were not considered in the calculation of the RMS error. Progression 
to the final surface RMS included holography system checkout, so does not 
represent the time required to set the VertexRSI antenna to its final surface. 



2 ) General Surface Stability: To estimate the accuracy and 
repeatability of the measurements, we produced difference 
maps between successive measurements throughout the mea- 
surement period. The RMS difference between consecutive 
maps is normally less than 10/im, typically 8/im. An example 



of a difference map is shown in Fig. 12 The map of measure- 
ment number 307 is shown on the left, while the right hand 
side shows the difference between map 307 and 308, made 
about one hour later. 

We have made many maps after the final setting of the 
surface while changing the orientation of the antenna with 
respect to the Sun. Maps taken over a few consecutive days 
were obtained at widely different temperatures, while also the 
wind conditions varied over time. The measured rms error 
during a 30 hour period in early May 2003 varied between 20 
and 23/im. During this period the wind was calm (< 5 m/s) 
and temperature changes of 15 C were encountered. A later 



4-may-2003-Vertex-0307 
T1S (unweighted) 23.45 /im 

ns (weighted) 19.75 /xm 

:vels: -100 to 100 by 20 fim 



diff 307 308 

rms (unweighted) 10.87 /xrr 

rms (weighted) 8.44 /xn- 

evels: -100 to 100 by 5 /irr 




Fig. 12. Example of the repeatability of the measurements. Horizontal and 
vertical axes units are meters with color scaling showing surface error in fim. 
The map on the right is the difference between the one at left and a map 
made one hour afterwards. The RMS of the difference maps is about 8/im, 
which is commensurate with the expected value due to noise and atmospheric 
fluctuations. 



5 day series in mid June gave rms errors from 22-26/im with 
temperature variation up to 20 C, wind velocities up to 10 m/s 
and periods of full sunshine. Most of this increase is believed 
to be due to the deteriorating atmospheric conditions at the 
VLA site during summer, when the humidity is significantly 
higher than normal. The much better results of 17/im obtained 
during the cold and dry winter period also point to a significant 
atmospheric component in the spring and summer results. 

However, some of the changes will be caused by tem- 
perature and wind. To increase the rms from 20 to 22/im, 
the "additional" component has a magnitude of 9/im rms. 
Such a contribution can be expected from the calculated 
values of Ajim each for wind and temperature for the panels, 
and 5/im for wind and 7/im for temperature for the BUS. 
These numbers are all within the specification. Actually, the 
measured differences are close to those expected from the 
estimated accuracy of the holography measurement and the 
measured rms differences in consecutive maps of about 8/im. 

B. AEC Antenna 

1 ) Overview: The apex structure of the AEC antenna does 
not enable us to mount the holography receiver inside the 
cylinder, as in the case of the VertexRSI antenna. Thus in this 
case the receiver was bolted to the flange on the "outside" 
of the apex-structure. Consequently, the feedhorn was brought 
to the required position by a piece of waveguide of about 
500 mm length. This caused significant attenuation in the 
received signal from the reflector to the mixer. Considering 
the available transmitter power, we concluded that this would 
not jeopardise our measurement accuracy significantly. 

The AEC antenna surface was set by the contractor with 
the aid of a Leica laser-tracker. The RMS of the surface was 
reported by the contractor to be 3 8/im. After this measurement 
a servo error caused the elevation structure to run onto the 
hard stops at high speed. The contractor decided to repeat 
the surface measurement and obtained an RMS of 60/im with 
some visible "astigmatism" in the surface. 

Our first holography map indicated an RMS of 55/im with 
a clearly visible astigmatism. We could identify the high and 
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low regions with those on the final AEC measurement. With 
two complete adjustments we surpassed the goal of 20/im. A 
third partial adjustment improved the surface RMS to about 
14/im. There is no indication of the "artefact" seen in the 
VertexRSI antenna. The results of the consecutive adjustments 



are summarised in Fig. 13 The last panel in this figure shows 
the final map after the repeated measurement and setting in 
January 2005. 





27-jan-2004 
52.6 i^m 



31 -jan-2004 
1 3.4 /^m 



03-feb-2004 
13.9 |im 



|l 

28-jan-2005 
17.5 /2m ^ 



Surface Error Map (/2m) 



-200 -100 100 200 
Error Distrib. (%) vs. Error (/xn 



Fig. 13. Sequence of surface error maps with intermediate panel setting. 
The surface contours are shown on the left side; the error distribution on the 
right. The white cross and the small white area represent the quadripod and 
a faulty panel and were not considered in the calculation of the RMS error. 
Note that the differences between the final three maps is within the overall 
uncertainty of the holography measurements. 

The adjustments were done with a tool provided by the 
contractor. It was similar to the one used by us on the 
VertexRSI antenna, but it was calibrated in "turns" rather than 
in micrometres. 

2) General Surface Stability: In Fig. 14 we show one 



of the final results and a difference map of this and the 
following measurement, made one hour later. The difference 
map indicates a repeatability of ~ 5/im RMS. There is no 
indication of the "ringing" in the outer region of the aperture, 
as was the case for the VertexRSI antenna. We ascribe this 
to the lower signal level due to the long piece of waveguide 
between feed and mixer. 



02-feb-2004-AEC-1 768 

rms (unweighted) 17.26 //m 

rms (weighted) 13.95 firu 

evels: -100 to 100 by 20 /im 



diff 1768 1729 

rms (unweighted) 6.17 jj^rr 

rms (weighted) 4.23 /xn- 

levels: -100 to 100 by 5 /irr 




Fig. 14. Example of the repeatability of the measurements. The map on 
the right is the difference between the one at left and a map made one hour 
afterwards. The RMS of the difference maps is about 5/im. 



Also here we made a series of 16 maps over a period of 
more than two days in early February 2004. Temperatures 
ranged from +2 to —10 C, while the wind was mostly calm 
with some periods of speeds up to 10 m/s. During one day 
there was full sunshine. The measured RMS error is very 
constant with a peak to peak variation of less than 2/im on an 
average of 14/im. The differences are fully consistent with the 
allowed errors under environmental changes and also of the 
same order as the measurement accuracy. We believe that the 
significantly better overall result is mainly due to the much 
drier and stabler atmosphere during these measurements as 
compared with the summer data from the VertexRSI antenna. 

VII. Conclusions 

We have successfully performed a holographic measurement 
and consecutive panel setting of the reflectors of the two 
ALMA prototype antennas to an accuracy of better than 20 
/im. Our estimated measurement accuracy is approximately 5 
/im. The data collection and analysis software packages are 
easy to use and provide quick results of the measurements, 
directly usable for a panel adjustment setting. We consider this 
system suitable for the routine setting of the ALMA production 
antennas to the goal of 20 /im accuracy in an acceptable time 
span. Modern survey equipment enables contractors to deliver 
reflectors with an accuracy of 50-60 /im without undue cost. 
Although the holography system can easily start with a much 
larger error, in the former case it is feasible to reach the 
specification with only one panel setting based on hologra- 
phy. We note that these measurements, being performed at 
one elevation angle only, do not provide information on the 
gravitationally induced deformation as function of elevation 
angle. 

In summary: 

1) The holography system has functioned according to 
specification and has enabled us to measure the surface 
of the antenna reflector with a repeatability of better than 
10/im. 



2) As shown in Figs. 11 and 13 we have set both antenna 
surfaces to and accuracy of 16-17 /im RMS. This will 
provide an aperture efficiency of about 65 percent of that 
of a perfect reflector at the highest observing frequency 
of 950 GHz. 
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3) The small differences in the surface maps obtained 
over several days of measurement are consistent with 
the measurement repeatability and at best marginally 
significant. If taken at face value, they indicate that 
the deformations of the reflector under varying wind 
and temperature influence are fully consistent with, and 
probably well within, the specification. This excellent 
behaviour over time is more important than the actual 
achieved surface setting. We stopped iteration of the 
settings after having achieved the goal of less than 20 
fim. 

4) Further information on the performance of the ALMA 
Prototype Antennas can be found in [13]. 



Appendix I 

Holography Map Parameter Equations and 
Calculations 



In this appendix we list the equations used to derive 



the holography measurement values listed in Table |III| The 
calculation leading to the power-related expressions of this 
appendix are detailed in [20]. 



/i = taper factor for signal feed 

= 1.13(6 + 2.5 dB taper) (28) 

fapo = apodization smoothing factor (29) 

fosr = map oversampling factor between rows (30) 

foss = map oversampling factor along a row (31) 

D = main antenna diameter (32) 

d = reference feed diameter (33) 

Oext = angular extent of map (assumed square) (34) 

Ob = primary beam size 

he 

vD 
61836.6/1 



v{GHz)D{m) ^^"^^^ 
sampling interval between rows 

A. 

fosr 

he 



(35) 



arcsec 



fosr^D 

61836.6/1 
fosMGHz)D{m) 
sampling interval along a scan 

foss 
Otsamp 

0.012^ arcsec 



(36) 



Nr. 



Lm 
P 

Pr 



6z 



B 

tint 



a 
R 
Az 



number of rows in map 

spatial resolution on dish 
D 



(38) 



flfa 



cm 



1717.7/1 /„po 
v{GHz)e^Adeg) 
map row scanning rate 
linear size of map 
Transmitter EIRP 

Reference feed power received 

7rd2 p 



(39) 

(40) 
(41) 
(42) 



4 47ri?2 
16 \R 

Main antenna power received on boresij 
7rZ)2 P 



4 

1 

16 



47ri?2 
R 



Ps{0) 



J ( TvaD ^ 

^1 y^T) 



2A 
Pr 



Ps {a)] kTsys 



sy ^ av 



16^2 fl Mo 

fos ^0 

Detector bandwidth 
Integration time 

■^row^row 
q2 



(43) 

;ht 

(44) 
(45) 
(46) 



(47) 
(48) 



fosr^Q 



1717.7 X 10'' fosr fiDjm) 
0{arcsec/sec)u{GHz)6'^{cm) 



hours (49) 



Scan angle, which ranges over ± (50) 
Distance between holography Tx and Rx(51) 
Reflector surface displacement accuracy (52) 



(37) 



For the ALMA holography system: 

1) tint = 36 msec, 

2) Mo = ^Ps{0)Pr = 4.167 X 10-^P, 

3) (Jo = (1.23 X 10-22w(P))^ 

4) Pr Term = (2.13 x 10-''^W{P)Y 

5) Average map noise for complex correlator {(lav) = 
(2.23 X 10-'^W{P))K 

)XlO~^ 

VP 



6) Sz = 
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Thus, if we want an error in the measurement of the surface 
shape of 5z = 5 fim, we need a transmitter with an EIRP of P 
= 7.3 /iW. The expected radiated power is in excess of 10// W, 
so there is a good margin. Noise will not be the limiting factor 
in the accuracy of the measurement. 

Appendix II 
Surface RMS Calculation Details 

The RMS in the holography plots is computed in the 
following way: 

• Unweighted: 




- The summation is on the N unmasked pixels 

- di is the normal surface displacement for pixel i 

• Weighted: 




- The summation is on the unmasked pixels 

- Ci is half of the path-length error due to the surface 
displacement for pixel i. This is directly related to 
the observed phase errors S(j)i by: 



47r 

= di cos ai (55) 

- Wi is the illumination amplitude of an ideal ALMA 
receiver at pixel i. This is currently specified to a 
12dB taper {wi is 0.251 at the edge of the dish). 
The function used is parabolic: 

= (1-10-0'^) (I)' (56) 
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